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ABSTRACT 

Context. Least-squares deconvolution (LSD) is a powerful method of extracting high-precision average line profiles from the stellar 
intensity and polarization spectra. This technique is widely used for detection, characterization, and detailed mapping of the temper- 

■ ature, magnetic, and chemical abundance structures on the surfaces of stars. 

Aims. Despite its common usage, the LSD method is poorly documented and has never been tested with realistic synthetic spectra. In 
Oh' this study we revisit the key assumptions of the LSD technique, clarify its numerical implementation, discuss possible improvements 

and give recommendations of how to make LSD results understandable and reproducible. We also address the problem of interpreta- 

■ tion of the moments and shapes of the LSD profiles in terms of physical parameters. 
Methods. We have developed an improved, multiprofile version of LSD (iLSD) and have extended the deconvolution procedure 
to linear polarization analysis taking into account anomalous Zeeman splitting of spectral lines. The iLSD method is applied to the 
theoretical Stokes parameter spectra computed for a wide wavelength interval containing all relevant spectral lines. We test various 
methods of interpreting the mean profiles, investigating how coarse approximations of the multiline technique translate into errors of 
the derived parameters. 

Results. We find that, generally, the Stokes parameter LSD profiles do not behave as a real spectral line with respect to the variation 
of magnetic field and elemental abundance. This problem is especially prominent for the Stokes / (intensity) variation with abundance 
and Stokes Q (linear polarization) variation with magnetic field. At the same time, the Stokes V (circular polarization) LSD spectra 
closely resemble the profile of a properly chosen synthetic line for the magnetic field strength up to 1 kG. The longitudinal field 
estimated from the Stokes V LSD profile is accurate to within 10% for the field strength below 5 kG and to within a few percent for 
the fields weaker than 1 kG. Our iLSD technique offers clear advantages over the standard LSD method in the individual analysis of 
different chemical elements. 

Conclusions. We conclude that the usual method of interpreting the LSD profiles by assuming that they are equivalent to a real 
spectral line gives satisfactory results only in a limited parameter range and thus should be applied with caution. A more trustworthy 
^s^j , approach is to abandon the single-line approximation of the average profiles and apply LSD consistently to observations and synthetic 

. spectra. 
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1 . Introduction Early applications of the multiline spectrum analysis focused 

on the detection of the weak magnetic field in cool stars us- 
Many modern echelle spectrographs allow one to record high- ing circular polarization observations. Direct detection of corn- 
resolution stellar spectra for the entire wavelength range span- plex and generally weak polarization signatures of the dynamo- 
ning from the near-UV to near-IR (3500-10000 A). These data generated mag netic fields in indi vidual spectral lines requires a 
contain an enormous amount of information about the physical S /N « 1000 dDonati et al.ll 19921). whic h is impossib le to reach 
conditions, dynamics, and chemistry of the stellar surface layers, for all but the brightest objects. ISemell (1 19891) and ISemel &~0 
However, often only a small number of spectral features is used ( 1996) suggested to average polarization measurements of many 



X 

& ' in astrophysical analyses that are based on these observations, spectral lines to detect and analyze ma gnetic fields with p olar- 

This is justified when one is interested in specific diagnostic lines ization observations of a modest quality. Donat fetaD d 19971) re- 

or in chemical elements with sparse spectra. On the other hand, fined this simple line addition procedure into the least-squares 

there are many processes that produce a very similar imprint on deconvolution (LSD) method. This technique assumes that all 

many spectral lines. For instance, the radial velocity shifts owing spectral lines have the same profile, scaled by a certain factor, 

to the orbital motion in a binary system, non-radial stellar pul- and that overlapping lines add up linearly. These assumptions 

sations, and cool spots on the surfaces of late-type active stars lead to a highly simplified description of the intensity and cir- 

affect all spectral lines in a similar way. Thus, one can employ cular polarization observations in terms of the convolution of a 

a multiline technique to extract common physical information known line mask with an unknown profile. With this approxima- 

from many spectral lines and determine an average profile with tion one can reconstruct an average line shape, the LSD profile, 

a high signal-to-noise (S /N) ratio. which is formally characterized by an extremely high S/N ratio. 

The LSD analysis of the circular polarization observa- 

Send offprint requests to: O. Kochukhov, tions has been very successful in detecting weak magnetic 

e-mail: oleg.kochukhov@fysast.uu.se fields in stars over the entire Hertzsprung-Russell diagram 
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(e.gJPetit et al.ll2008HLignieres et al.ll2009tlAuriere et all2 009: 



Alecian et alJl2008t iDonati et al]l2008t iMorin et al.ll2008l) . The 
LSD was also ext ended to the Stokes Q and U spectra by 
Wade et a"D (|2000b) with the aim to characterize weak linear po- 
larization signatures in magnetic Ap stars. Outside the research 
area of stellar activity and magnetism LSD was applied to dis- 
entangle composite spectra of multiple stars and me asure their 
orbital radial velocity variation (lHareter et al.l 12008). A modi- 
fication of the LSD method was also employed in the search 
for stellar differential ro tation through the line profile analysis 
dReiners & Schmittll2003l) . 

Many studies went further than using LSD as a simple 
magnetic field and starspot detection tool, attempting to in- 
terpret the LSD profiles and extract various physical infor- 
mation. These applications range from a relatively simple de- 
terminatio n of the long i tudin al magnetic field from the LSD 
Stokes V (lAuriere et al.l 120071) and of the net linear p olariza- 
tion from the LSD Stokes QU data dWade et al.ll2000al) to nu- 
merous sophisticated Doppler and Zeeman-Doppler imaging 
studies. In the latter investigations the time series of the LSD 
profiles were interpreted under different physical assumptions, 
with regularized in version tech niques to produce magnetic (e.g., 
iDonati etai] 120031). chemica l dFolsom et al.ll2008l) . and bright- 
ness (e.g. JJeffers et al.ll2007l) maps of stellar surface. All these 
studies implicitly assumed that the LSD profile behaves similar 
to a real, isolated spectral line with average parameters. 

The proliferation of the studies that use LSD outside its origi- 
nal validity range and assign certain meaning to the LSD profiles 
without theoretical justification or numerical tests casts doubt 
on some of the results obtained with this method. Few p apers 
attemp ted to clarify or improve the multiline technique. iBarnesI 
(2004) studied mathematical aspects of the spectral deconvolu- 
tion procedure. Sennhauser et al. (2009) attempted to improve 
the LSD description of the intensity spectrum by accou nting 
for nonlinearity in blende d profiles. ISemel et' al. (2009) and 
Rami rez Velez et al] (1201 Ol) tested alternative multiline Zeeman 
signature methods (similar but not identical to LSD). 

None of the previous studies came close to characterizing the 
LSD technique and testing the reliability of the LSD profile in- 
terpretation for a realistic situation of blended intensity and po- 
larization profiles of thousands of spectral lines. This task was 
beyond the capabilities of a physically complete spectrum syn- 
thesis for spotted magnetic stars, which was typically limited 
to a few dozen spectral lines in s everal short wavelength inter- 
vals (e.g., Kochukhov et al. 2004). Recent progress in numerical 
techniques of the polarized spectrum synthesis and continuous 
expansion of computer resources has allowed us to overcome 
these limitations. In this paper we present the first comprehen- 
sive, direct numerical tests of the LSD technique, complemented 
with an improvement and clarification of the LSD method it- 
self. We compute four Stokes parameter stellar spectra in a wide 
wavelength region, including all relevant metal absorption lines. 
We then apply LSD to these synthetic spectra and compare quan- 
tities inferred from the LSD profiles with the input parameters 
adopted for polarized spectrum synthesis. This work allows us 
to reach general conclusions about the validity of different ap- 
plications of LSD in stellar physics. 

Our paper is structured as follows. In Sect. [2] we revisit the 
general principle and basic assumptions of the LSD technique. 
Then we describe numerical implementation of the deconvolu- 
tion procedure in our LSD code and consider possible improve- 
ments of the least-squares deconvolution technique. We also ad- 
dress the question of normalization of the LSD profiles and out- 
line several methods of their interpretation. Sectionals devoted 



to numerical tests of the LSD technique with the help of theoret- 
ical polarized spectrum synthesis. We start with a brief descrip- 
tion of the magnetic spectrum synthesis code that we use to sim- 
ulate realistic spectra in four Stokes parameters. We examine the 
behavior of the LSD profile shapes and of the integral quantities 
that can be inferred from the average profiles in different Stokes 
parameters (longitudinal magnetic field, net linear polarization, 
equivalent width). The paper concludes with the summary of the 
main results of our investigation in Sect. [4] 



2. Least-squares deconvolution 

2.1. General principle 

We consider a situation where each line in the intensity or polar- 
ization spectrum of a star can be approximately described by the 
same line shape, scaled by a certain factor. If we further assume 
that overlapping line profiles add up linearly, the entire spectrum 
Y(v) can be described as a sum of scaled and shifted identical 
profiles Z(v) 



(1) 



where Y is a model residual intensity, 1 - I/I c , or normalized po- 
larization spectrum, e.g. V/I c . The position in the velocity space 
V,- = cAAi/Ai is associated with a wavelength shift AAj from the 
central wavelength of the z'-th line. The relative contribution of 
each spectral l ine is given by the w eight w, . 

Following Don ati et"a"D ( 19971) . Eq. (Q} can be expressed as a 
convolution of the line pattern function 



Af(v) = Wi5(v - vj) 



and the mean profile Z(v) 



Y = M*Z, 



or, equivalently, as a matrix multiplication 
Y = M ■ Z, 



(2) 



(3) 



(4) 



where Y is an n-element model spectrum vector, Z is an m- 
element common profile spanning a certain range in velocity, 
and M is a n x m line pattern matrix containing information on 
the line positions and relative strengths. In typical applications 
m might be on the order of 10 2 or 10 3 , while n might be on the 
order of 10 5 or 10 6 . 

The scope of the LSD technique is to solve an inverse prob- 
lem corresponding to Eq. (O or |@}. Namely, to estimate the 
mean line profile Z for a given line pattern information M, ob- 
served spectrum Y° and associated statistical uncertainty cr. In 
general, this approach is justified as long as the average pro- 
file can be interpreted meaningfully in terms of useful physical 
quantities, and its uncertainty <x z is much smaller than the error 
of the input spectrum, cr z <k cr. 

An average line profile can be also estimated by the direct 
superposing, weighting, and averaging a large number of short 
pieces of spectrum. This approach would yield results identi- 
cal to LSD if all lines were free of blends. In reality, absorption 
lines in stellar spectra are rarely isolated. For this reason direct 
averaging recovers a profile that is systematically distorted by 
overlapping lines, whereas LSD is able to provide a clean mean 
profile. 
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Note that in its application to the Stokes / spectrum the LSD 
method is similar to the broadening fu nction (BF) analysis tech- 
nique developed by iRucinskil (|1992|) . The essential difference 
is that the BF technique starts from an unbroadened observed 
or synthetic template spectrum rather than a superposition of 6- 
functions employed by LSD. The requirement of using a realis- 
tic input template imposes certain limitations in practical appli- 
cations of the BF method and does not allow a straightforward 
extension of this technique to the multiline polarization analy- 
sis. On the other hand, the BF method provides a more precise 
approximation of the intensity spectrum because it does not rely 
on the assumption of a linear addition of blended lines in the 
template spectrum. 

2.2. Physical assumptions 

The simplified description of the stellar spectrum adopted by 
LSD seems to be very coarse and is obviously invalid in some 
extreme cases. For instance, there is very little similarity of the 
spectral line profiles when the magnetic field strength is so large 
that individual Zeeman splitting patterns of spectral lines are 
resolved. But in some other cases, such as unsaturated lines in 
weak magnetic field, the LSD assumptions are not so bad. The 
goal of this section is to discuss which physical approxima- 
tions correspond to the mathematical expressions employed in 
the LSD method and to qualitatively estimate the validity range 
of these assumptions. 

2.2.1. Intensity 

Let us start from the assumption of the line profile similarity in 
non-magnetic case. It is known that if the lines are unsaturated, 
the local residual intensity profile emerging from a stellar atmo- 
sphere is proportional to the line absorption coefficient m at a 
certain depth 

1 - X' e (A)/X' c cc K[ (A) (5) 

and, therefore, can be represented as an average local profile Zj^ 
scaled in depth by the corresponding central line depression d\oc 
and in width by the wavelength A 

l-X'^/X'^d^Z'^v). (6) 

Thus, for LSD to be applicable to the local Stokes / spectrum 
X 1 , spectral lines must all have a similar shape of the absorp- 
tion coefficient K( and be relatively weak. The former condition 
suggests that the multiline technique should avoid spectral fea- 
tures with an unusual line shape, such as the H i and He n lines in 
hot stars or metal lines exhibiting strong damping wings in cool 
stars. 

A simple assessment of the validity of the line shape self- 
similarity condition can be obtained by considering realistic the- 
oretical profiles of typical spectral lines of a different strength. 
With the spectrum synthesis technique described in more detail 
below (Sect. 13.11 ). we calculated local profiles of fictitious Fe i 
line with representative parameters and oscillator strength ad- 
justed to yield central depths from 5% to 90% of the continuum 
intensity. The discrepancy between the real line profile and a 
scaled version of the weakest line is illustrated in Fig. [T]for an 
infinite resolving power and for R = 60 000. Keeping in mind 
the typical observational uncertainty of ~ 10~ 2 of the modern 
echelle spectra of bright stars, one can conclude that an identi- 
cal line shape assumption of the LSD of Stokes / is appropriate 
only for lines weaker than ^40%. These results were obtained 




Line depth 

Fig. 1. Discrepancy between scaled profile of the weak Fe i spec- 
tral line and detailed spectrum synthesis calculation for the line 
of given central depth. The discrepancy is characterized with the 
standard deviation (solid line) and maximum difference (dashed 
line). Results are shown for an unbroadened spectrum (circles) 
and for the resolving power R = 60 000 (squares). 
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Line depth 



Fig. 2. Same as Fig.[TJfor the discrepancy between linearly added 
profiles of two identical Fe i lines of a given central depth and 
detailed spectrum synthesis calculation. 



with a model atmosphere typical of a mid- A star (T e g = 9000 K, 
log g = 4.0), but a calculation with a solar model atmosphere 
yielded very similar results. 

An analogous calculation can give us an idea about the ac- 
curacy of the linear addition of the blend components. This ap- 
proximation is evidently valid for the spectral lines described 
by Eq. (O, but it becomes increasingly unrealistic for stronger 
lines. Figure|2]illustrates the discrepancy between detailed spec- 
trum synthesis of the blend consisting of two identical Fe i lines 
with a given strength and simple arithmetic addition of the line 
depressions corresponding to the individual blend components. 
This figure shows that the accuracy of the linear line-adding re- 
mains at the level of < 10~ 2 if the residual intensity of the blend 
components is S= 20%. 
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Until now we were examining the behavior of the local in- 
tensity profiles which, of course, are not directly observable in 
real stars. In reality, the LSD method approximates the disk- 
integrated stellar spectra. Integration over all points M on the 
visible stellar disk takes into account the distribution of the con- 
tinuum brightness biw and rotational Doppler shift vm 



1 - I(v)/I c 



I 



b M d\, 



0C Z[ 0C (M,v 



v M )dS * dZ'(v). 



(7) 



Here d is the line depth in the disk-integrated spectrum and the 
second equality is approximate because we have assumed that 
the center-to-limb variation of the normalized local line profiles 
of all lines can be described by the same linear law, independent 
of the wavelength and atomic line parameters. This assumption 
is not a particularly restrictive one for the majority of normal 
stars. However, one can find examples, such as the spectra of 
rapidly rotating early-type stars, where the line-to-line difference 
in the center-to-limb variation is emphasized by an inhomoge- 
neous temperature distri bution and can b e readily detected in the 
stellar flux spectrum (Takeda et al. 200H). 

2.2.2. Polarization 

The characteristic splitting of the Zeeman components of a spec- 
tral line is 



AA P _ 



A 2 c e Q B 



(8) 



where A c is the wavelength of the unperturbed line, B is the local 
field strength, and the other variables have their usual meaning. 
When the magnetic splitting is much smaller than the intrinsic 
width of the spectral lines determined by the thermal Doppler 
broadening A/Id, one can derive approximate expressions relat- 
ing intensity (Stokes I) with the Stokes parameters characteriz- 
ing circular (Stokes V) and linear (Stoke s QU) polarization (e.g., 
lLandi Degl'Innocenti & Landolfill2004l) . 

In the weak field limit, formally defined by the condition 
AA B <s AAo, any spectral line obeys 



..... _A/1b dl 
V(A)CC8 AA^dA- 



(9) 



The effective Lande factor g entering this expression is given by 

g= +#2) + jigi -gz)r, (10) 
where 



r=J l {h + \)-J 2 {J 2 +\), 



(11) 



This relation between Stokes V and the first derivative of Stokes 
/ is valid for lines of any strength. The corresponding weak-field 
relation for the linear polarization, 



Q,UocG\ 



'AAk\ 2 d 2 I 



AA D dA 2 



with 



1 



G = g 2 -—( gl -g 2 ) 2 (16s-7r 2 -4) 
and 

s = /i(/i + 1) + J 2 (J 2 + 1), 



is derived using a set of more restri ctive assumptions (see 
lLandi DeglTnnocenti & Landolfil 120041) . essentially equivalent 
to the weak line approximation of Sect. 12.2.11 

Following the same arguments as given above for the trans- 
formation from the local to disk-integrated intensity spectra and 
taking into account that AAg/AAr, ~ A c , we conclude that the 
disk-integrated Stokes V profiles can be represented by a com- 
mon profile scaled by a factor wy = dA c g: 



V{v)jl c * w v Z v (v). 



(15) 



The linear polarization line shapes are approximated with 
Q{v)lh * w Q Z Q (v), (16) 

wh ere the linea r polariz ation line weight wq = dA 2 G. 

IWade et al.l d2000bl) established a similar relation for the 
scaling of the linear polarization profiles, but they used g 2 in- 
stead of G in the line weight wq. This is equivalent to assuming 
that all spectral lines have a simple triplet splitting pattern. In 
general, G < g 2 for the anomalous Zeeman splitting, and some- 
times G < 0. In Sect. [3] we demonstrate that our generalization 
of the LSD linear polarization profile scaling factor yields a bet- 
ter approximation of the Stokes parameter spectra obtained with 
detailed spectrum synthesis calculations. 

We emphasize that the self-similarity of the polarization pro- 
files asserted by Eqs. (15[ and ( TT~6b is derived assuming both the 
weak field regime for the magnetic splitting and the weak line 
approximation for the line formation (implicitly for Stokes V and 
explicitly for Stokes QU). The validity of the second assump- 
tion was assessed in Sect. 12.2.11 On the other hand, the weak 
field approximation does not have a well-defined validity range. 
Most authors agree that it is fairly accurate only for the mag- 
netic field strengths below ~ 1 kG. However, this approximation 
is also w idely used for the LSD analysis of multi-kG fields of 
Ap stars dWade et al.ll2000al) . 

We suggest that it makes little practical sense to assign a gen- 
eral "validity range" to the LSD polarization profiles and quan- 
tities derived from them. Instead we aim to quantify the LSD 
profile behavior numerically and assess the robustness of their 
interpretation in the context of specific applications, taking into 
consideration the quality of observational data. 

2.3. Numerical implementation 

We follow iDonati et alJ (1 19971) in representing the LSD model 
spectrum as a matrix multiplication described by Eq. ©. We 
use linear interpolation to project the mean profile Z(v) onto an 
arbitrary wavelength grid of observations. Then a spectral line 
with a weight Wi and central wavelength Ai contributes along a 
bidiagonal of the line mask matrix M 



with 



Wl(.Vj+\ - Vi)/(Vj + l - Vj) 

wiivj - Vj)/(v j+ i - vj) 



(12) vi = c(Ai - Ai)/A/ and vj < v-, < v j+i . 



(17) 



(18) 



Here index i corresponds to the wavelength space of observa- 
tions, whereas index / runs over all velocity points of the LSD 
profile. 

(13) The problem of finding the best-fitting mean profile Z(v) for 
a given line mask M and observations Y° is equivalent to mini- 
mizing the x 2 function 



(14) x 2 - (Y° - M ■ Z) T ■ S 2 ■ (Y° - M ■ Z) — > min, 



(19) 
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where S is the square diagonal matrix containing inverse vari- 
ance, S a = 1 /o~i, for individual pixels in the observed spectrum. 

Since the model describing observations is linear, the least- 
squares solution of the x 1 problem formulated in Eq. ( fT9l l is 
given by 

Z = (M r ■ S 2 • Mr 1 ■ M r • S 2 ■ Y°. (20) 

The M r ■ S 2 ■ Y° term in Eq. ( l20b represents a weighted cross- 
correlation between the line mask and the observed spectrum. 
The left-hand part is the inverse of autocorrelation matrix M r • 
S 2 • M, which effectively deconvolves the cross-correlation vec- 
tor and provides formal uncertainty estimates of the LSD profile 
through the diagonal elements of (M r • S 2 • M)" 1 . 

We have implemented the numerical technique described 
above, along with the improvements related to the linear po- 
larization treatment (Sect. 12.2.21 ). multiple mean profile analysis 
(Sect. 12.4.11 ), and regularized deconvolution (Sect. 12.4.21 ). in an 
improved least-squares deconvolution (iLSD) Fortran code. On 
input it reads intensity or polarization observations and the cor- 
responding error bars. Several sets of line weights are defined ex- 
plicitly in another input file. Then, the code fills the line pattern 
matrix M for a prescribed mean profile velocity grid and evalu- 
ates Z(v) using Eq. d20l >. The inverse of the autocorr elation ma- 
trix is calculated with the help of LU-decomposition dPress et alj 
Il992h . The error bars of the resulting LSD profile are scaled to 
ensure that the reduced chi-square is equal to 1. In our imple- 
mentation of the LSD technique we consider only those pixels 
of the observed spectrum that contain the contribution of at least 
one spectral line. 

As discussed in Sect. 12.2.11 the coarse LSD model spec- 
trum calculated with Eq. (0]) significantly overestimates the in- 
tensity of strong ble nds, sometimes leading to unphysical results. 
iDonati et all (1 997) mention that they ensure that "the sum of 
normalized central depths over nearby lines never exceeds 1". 
However, it is not clear how this constraint can be implemented 
within the linear LSD approach. G. Wade (private communica- 
tion) informs us that this modification is not used in the practi- 
cal applications of Donati's LSD code. Thus, we also solve the 
least-squares problem defined by Eq. ( fT9l ) without imposing any 
additional constraints. 

2.4. Improvements of deconvolution technique 
2.4.1. Multiprofile LSD 

A straightforward generalization of the single-profile LSD tech- 
nique can be obtained by postulating that each line in the stellar 
spectrum is represented by a superposition of N different scaled 
mean profiles: 

N 

Y = Y J M k -Z k =M'-Z'. (21) 
k=l 

In that case one can still use the matrix formulation of the LSD 
technique, but with a composite line pattern matrix M' with the 
dimensions n x (m ■ N). The matrix M' is formed by concate- 
nating matrices that have the same structure, but a different 
magnitude of the non-zero elements corresponding to N differ- 
ent sets of line weights. All N average profiles, concatenated in 
the composite profile Z', are reconstructed simultaneously with 
the matrix inversion technique described in Sect. 12.31 

There are several advantages of using multiprofile LSD. One 
can, to some extent, account for the line shape difference be- 
tween weak and strong lines or recover several mean profiles 
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Fig. 3. Top panel: Comparison between the usual LSD profile, 
the scaled cross-correlation function and the LSD profile recon- 
structed using regularization. These profiles are obtained from 
a simulated S/N=15 spectrum, containing 500 randomly dis- 
tributed scaled lines with a characteristic Stokes V S-shaped sig- 
nature. Profiles are shifted vertically for display purpose. The 
indicated signal-to-noise ratio is inferred from the pixel scatter 
in the regions on both sides of the profile. Bottom panel: The 
difference between LSD profiles and the corresponding CCF. 

if the stellar spectrum contains two or more systems of lines 
formed under sufficiently different physical conditions (e.g., a 
composite spectrum of a binary star). 

In the studies of magnetic chemically peculiar stars, which 
often exhibit diverse distributions of chemical ele ments, one 
has to c onsider LSD p rofiles of individual elements (Wade et al. 
I2000bt lFolsometal.[ |2"008). Using conventional single-profile 
LSD for a given element often introduces significant distortions 
in the wings of the corresponding mean profile owing to ne- 
glected blending by lines of other elements. In contrast, with the 
multiprofile least-squares deconvolution we are able to account 
for the background lines blending the element of interest. We 
have implemented multiprofile LSD in our iLSD code and apply 
it extensively in the numerical tests presented in Sect. [3] 

2.4.2. Regularized deconvolution 

With the multiline technique described by Eq. d20l one derives 
the mean profile in two steps. First, a cross-correlation function 
(CCF) is obtained and, second, it is multiplied by an inverse of 
the autocorrelation matrix to remove the spurious pattern aris- 
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ing from possible regularities in the line distribution. However, 
the latter operation suffers from a well-known problem of noise 
amplification, which is common to all applications of direct de- 
convolution. In the present case the second step systematically 
increases the noise in the mean profile compared to CCF inde- 
pendently of the numerical technique used to invert the auto- 
correlation matrix. This is illustrated in Fig. [3] which shows the 
reconstruction of the mean profile with a typical Stokes V signa- 
ture from 500 superimposed, scaled copies of this profile that are 
randomly distributed in wavelength. Random noise was added to 
simulate observational data with a signal-to-noise ratio of 15. 

A simple modification of the LSD procedure allows circum- 
venting the profile quality degradation associated with decon- 
volution. It also allows us to decouple the mean profile noise 
characteristics from the step size of the velocity grid on which 
the LSD profile is reconstructed. We modify Eq. fl% to require 



X + AT? — > min, 



(22) 



where "R is the regularization function, penalizing large pixel- 
to-pixel variations of the LSD profile, and A is the correspond- 
ing regularization parameter. The first-order Tikhonov regular- 
ization, 



* = 7 Tj i (Zi ~ z '- l)2 + (Z/ " Z/+l)2 ] ■ 



(23) 



which minimizes the first line profile derivatives, is a convenient 
choice of regularization for this problem. It can be implemented 
in the LSD matrix algorithm by modifying Eq. ( f20b as follows 



Z = (M T ■ S 2 ■ M + ARr 1 ■ M T ■ S 2 ■ Y°, 

where the tri-diagonal m x m matrix R is given by 



(24) 



R 



1 -1 
-1 2 



-1 2 
-1 



(25) 



Figure [3] demonstrates that with an appropriate choice of 
the regularization parameter one can deconvolve the cross- 
correlation function without degrading its noise properties. In 
this particular example the regularized deconvolution yields 
roughly 50% gain in S/N compared to the standard LSD tech- 
nique. 

The least-squares deconvolution constrained by the 
Tikhonov regularization of the mean profile is implemented in 
our code, but is not applied in the numerical experiments of 
Sect.[3]because those are conducted for the noise-free spectra. 

2.4.3. Nonlinear deconvolution 

Several studies suggested modifications of the LSD method 
that would overcome its most important limitation: the assump- 
tion of linear addition of the overlapping lines, which leads 
to a poor description of the Stokes / spectrum. The simplest 
way to reduce the discrepancy between the coarse LSD de- 
scription and real intensity spectrum is to adjust not only the 
mean profile, but also th e line strengths themselves, as done by 
Reiners & Schmitt (2003) in their "physical least-squares decon- 
volution". This enhancement produces a superficially better fit 
to Stokes /. However, the mathematical formulation and numer- 
ical implementation of this method is somewhat ambiguous be- 
cause one cannot determine the amplitude of the mean profile 



and line strengths in the LSD mask at the same time. One can 
also encounter uniqueness problems when fitting intensities for 
a group of very close lines. In addition, this technique cannot 
be easily extended to polarization analysis. In general, it yields 
the adjusted line strengths that are valid for Stokes / but may be 
inappropriate for other Stokes parameters. 

Recognizing that the main li mitation of the classical LSD is 
its primitive treatment of blends. ISennhauser et al.l d2009t) devel- 
oped "nonlinear deconvolution with deblending" (NDD) tech- 
nique in which the overlappin g lines ar e added nonlinearly, using 
a formula derived from the Minnaert d 19351) approximation of 
the lin e shapes in the solar spectrum. Although lSennhauser et al.l 
(2009) argued that NDD is superior to LSD in the description 
of Stokes /, in their method they failed to distinguish two dif- 
ferent types of blending present in the stellar spectrum. Very 
close spectral lines with overlapping absorption coefficients in- 
deed add up nonlinearly. But the residual intensities of separated 
lines whose profiles overlap because of rotational, macroturbu- 
lent or instrumental broadening are simply added together, just 
like the standard LSD assumes. In Fig. [4] we illustrate the differ- 
ence between NDD and LSD description of blends in these two 
cases. Evidently the NDD method does not always yield a better 
approximation of the blended lines and in many cases appears to 
be inferior to LSD. 

We agree that a nonlinear treatment of overlapping lines is 
a promising direction for improving LSD. However, one has to 
somehow distinguish the linear and nonlinear blending. One can 
envisage a successor to LSD that uses a method similar to NDD 
for the calculation of an unbroadened tem plate spectrum a nd 
then uses the broadening function approach (RucinskU l 19921) to 
account for various broadening effects that can be treated as con- 
volution. We postpone the development and implementation of 
this two-stage multiline technique to future papers. 

2.5. Normalization of the LSD profiles 

A mathematical description of the stellar spectra by the LSD 
technique does not allow one to determine the mean profile 
amplitude independently of the scaling of weights in the line 
mask M. Observations in each Stokes parameter constrain the 
product of the corresponding line weight and the LSD profile. 
Consequently, the mean profile amplitude depends on the arbi- 
trary normalization of the line weights. This normalization deter- 
mines parameters that should be assigned to the resulting mean 
profile. For example, the general form of the LSD line weights 



dAg 



doA g 
for Stokes V and 



w QU 



dA 2 G 
doA^Go 



(26) 



(27) 



for linear polarization. 

Normalization by the parameters fo = do, Ao,g and Go can 
be equivalently applied to the input LSD weights or to the re- 
sulting mean profi les. According to Sect. 4.1 of their paper, 
iDonati et al.l d 19971) have arbitrarily chosen Aq = 500 nm, do — 1 
and g Q = 1. In other words, their reconstructed LSD profiles 
should correspond to a line with both the central intensity and 
the effective Lande factor set to unity. Most other studies em- 
ploying the LSD technique neglected to explicitly describe the 
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Fig. 4. Comparison of the LSD (dashed line) and NDD (dash-dotted line) description of blended lines with the detailed radiative 
transfer (RT) calculation (solid line), a) Two identical unresolved lines with overlapping absorption coefficients, b) Two resolved 
lines that overlap because of the rotational broadening corresponding to v e sin i = 10 km s . c) Same as b), but the overlap is caused 
by the Gaussian instrumental broadening with FHWM = 0.17 A. 



adopted normalization of the LSD weights, thus obscuring the 
interpretation of their published mean profiles. 

On the other hand, many LSD-based investigations quote the 
mean line parameters, usually computed taking into account the 
statistical uncertainty of the observed spectrum: 



(/> = 



1,1 filer] 



(28) 



where // is the parameter that is averaged and o~i is an esti- 
mate of the observ ational uncertainty at the location of line I. 
Occasionally (e.g. jBarnes et al.|[20 00). the averaging takes into 
account the line weight itself 



<f)w 



(29) 



Compared to Eq. d28l i, this method yields systematically dif- 
ferent average parameters because fi and w/ are correlated. For 
example, for a sample of 8864 lines used for the numerical 
experiments in Sect. [3] of our paper we find (g) = 1.21 and 
(g)w - 1.41 for the weight wy computed with Eq. (f26b . Most 
LSD studies do not comment on the averaging method they use. 

Note that no matter which expression is adopted for com- 
puting the average line parameters, these quantities characterize 
only certain statistical properties of a set of lines included in 
the LSD mask, but shoul d not be assigned to the LSD profiles 
themselves as implied bv lDonati et al. d 1997b and many subse- 
quent studies. The parameters corresponding to the LSD profiles 
are determined uniquely by the arbitrary normalization constants 
/o, which are not necessarily the same as the average parameters 
computed with Eqs. d28l > and ( l29l i. unless a special effort is made 
to ensure fo = (/) or fo = (/)„„ as done by some LSD studies (V. 
Petit, private communication). Of course, the latter approach be- 
comes advantageous when one uses a spectrum synthesis code, 
trying to interpret the LSD profiles as Stokes parameters of a real 
spectral line (see Sect. 13.5b . 

The confusion between the average line parameters and the 
LSD weight normalization persists through many papers and has 
possibly led to erroneous interpretations of the mean profiles in 
some of these publications. We recommend that studies employ- 
ing the LSD method should follow simple guidelines in order 



to better document application of this technique and make their 
results transparent and reproducible: 

- provide values of the parameters do, Aq, g Q , and Go (if linear 
polarization is involved) adopted for normalizing the LSD 
weights or renormalizing the LSD profiles; 

- clarify which expression is used to compute the average line 
parameters; 

- most importantly, explain which meaning is assigned to the 
average quantities and state whether or not the presented 
LSD profiles are rescaled to ensure fo = (/). 

2.6. Interpretation of the LSD profiles 

The average intensity and polarization profiles provided by the 
LSD technique can be interpreted with various assumptions and 
levels of sophistication. In addition to merely using LSD as a 
qualitative magnetic field detection tool dDonati et al Jl 19971) . we 
can distinguish the following main directions of utilizing them 
in the stellar physics studies: 

1 . One can infer characteristics of a star and its magnetic field 
through a direct comparison of the LSD profiles obtained 
self-consistently (i.e. using the same line lists and weights) 
from observations and from the detailed theoretical spectrum 
synthesis of the entire stellar flux spectrum in several Stokes 
parameters. Then the LSD profiles are used only for the inter- 
comparison of theory and observations but are not assigned 
any particular physical meaning. Obviously, this is the most 
robust method of treating the output of the least-squares de- 
convolution, but it is still far beyond the capabilities of the 
current computational resources and thus hasn't been applied 
in any real study. 

2. Taking into account that the LSD procedure is linear, we 
can use the theoretical spectrum synthesis of the entire stel- 
lar intensity spectrum to derive local LSD profiles for a 
range of temperatures, chemical abundances, magnetic fields 
strengths, and orientations. Interpolating within this theoret- 
ical local LSD line profile table we can obtain the model 
LSD flux profiles in one or several Stokes parameters and 
compare them with the corresponding LSD spectra inferred 
from observations. 
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This analysis technique has not been implemented yet. 
However, some studies employed a simpler version of this 
method where the LSD profiles obtained from the observed 
flux spectra of sharp-lined benchmark stars of different T e ff 
are assumed to represent local profiles corresponding to dif- 
ferent regi ons on the surface o f spotted rapidly rotating ac- 
tive stars (Bar nes et al.l 12000). This approach neglects the 
variation of line profiles with the limb angle and typically 
relies on a linear interpolation between only two template 
stars. An extension to the circular polarization LSD analysis 
is achieved within the framework of the weak field approxi- 
mation, adopting derivatives of the LSD Stokes / profiles of 
sharp-lined stars as the local circ ular polarization spectrum 
dDonati & Collier Cameronll 19971) . 

3. Assuming that the LSD profiles behave as a real spectral line 
with average parameters, one can interpret the LSD spectra 
with standard theoretical tools. This approach is frequently 
used in recent Doppler imaging studies of the stellar sur- 
face structure and magnetic fields (e .g., iDonati et"ai]l2008t 
iMorin et ai]|2008t iFoisom et al.ll2008l) . 

4. Finally, we can limit ourselves to the analysis of moments 
of LSD profiles that are assumed to behave similar to the 
corresponding moments of real spectral lines. In particular, 

- The zeroth-order moment of Stokes /, equivalent width 



W v = do 



/Z'dv, 



(30) 



can be used to measure metallicity or abundances of in- 
dividual elements. Note that in the expression for W v 
the integral of Z 7 is scaled by the central line intensity 
c/() adopted for the normalization of the Stokes / LSD 
weights. 

- The first-order moment of Stokes V provides a measure 
of the mean longitudinal magnetic field 

£ f vZ v dv 

(B x ) = -7.145 xlO 6 J _ - , (31) 
'Wo J z dv 

where the magnetic field is measured in Gauss, the wave- 
length in A, and velocity in kms -1 . The wavelength Aq 
and effective Lande factor g appearing in this expression 
represent the arbitrary quantities adopted for the normal- 
ization of the Stokes V LSD weights (see Eq. (Got), but 
not the mean lin e parameters, as often erroneously sta ted 
in LSD papers dDonati et al.ll 19971: IWade et alJl2000ah . 

- The normalized equivalent width of the LSD Stokes Q 
and U LSD profiles 

( WQU )fzOd V (w QU )fZ u dv 

Pq = pr— and P v = — , (32) 

j Z'dv J Z'dv 

where (wqu) is the average Stokes QU L SD weight, 

charac terizes the net linear polarization (Wade et al. 

2000a) and has a similar information content as the 

broad-band linear polarizat ion measured with the help of 

photopolarimetric devices dLerov|[l995l) . 

It is beyond the scope of this paper to explore all these meth- 
ods of interpreting the LSD profiles. The numerical experiments 
presented below are mainly focused on testing the third approach 
and several aspects of the fourth approach. 



3. Numerical tests of the LSD technique 

3.1. Polarized spectrum synthesis 

We use the magnetic spectrum synthesis code synmast to calcu- 
late four Stokes parameter stellar spectra covering a large part of 



the optical wavelength range, synmast is a fur ther devel o pmen t 
of t he synthmag code, described in detail by Piskunov (1999) 
and lKochukhovl d2007l) . We refer the readers to these papers for 
an in-depth discussion of the numerical techniques and input 
physics implemented in the code. In summary, we solve numer- 
ically the polarized radiative transfer equation for a given stel- 
lar model atmosphere and homogeneous or vertically stratified 
chemical abundances assuming local thermodynamical equilib- 
rium. The code includes an upd ated version of the molecular 
equilibrium solver developed bv IValenti et al.l ([T998), allowing 
us to treat Zeeman and partial Paschen-Back splitting in both 
atomic and molecular absorption lines. 

Previous studies extensively used synthmag calculations of 
individual spectral lines or blends in narrow wavelength inter- 
vals to investigate the magnetic field, chemical composition, 
and stratification in A and B chemically pecul iar stars (e.g., 
lKochukhovll2003L 120061: Ikvabchikova et al^OOol) and to model 
magnet ic fields of T Tauri stars djo hns-Krull 2007) and active M 
dwarfs dJohns-Krull & Valentilll996l:lKochukhov et al.ll2009h . In 
the present study we have extended the capabilities of the code, 
adapting it for massive polarized spectrum synthesis of thou- 
sands of spectral lines. 

The output of synmast is a set of local Stokes IQUV spectra 
for a given magnetic field vectoiQ and the angle of view from 
the normal to the surface (yU = cos 9). These spectra can be com- 
bined with appropriate weights and Doppler shifts to simulate 
stellar flu x spectra for a sim ple homogeneous magnetic field dis- 
tribution (Kochukhov 20 071) or for an arbitrary pre scribed map of 
surface structures dPiskunov & Kochukhovl20()2l) . However, be- 
cause each of these flux spectra represents a linear superposition 
of the local intensity spectra computed for different magnetic 
field parameters and the LSD procedure behaves as a linear op- 
erator, it is sufficient to limit the investigation to the local, disk- 
center Stokes parameters that are calculated for a sufficiently 
broad range of the field strengths and inclinations. This analy- 
sis should provide a basic but still rather general assessment of 
the performance of the LSD technique without digressing to the 
analysis of particular surface magnetic field configurations. 

For the magnetic and chemical abundance LS D experiments 
presented in this section we considered an atlas9 (lKuruczll 993) 
model atmosphere of the star with T e ff = 9000 K and log g = 4.0. 
We adopted solar abundances for all elements except Cr and Fe. 
The concentration of these two species was enhanced by 2 and 
1 dex, respectively, which is typical of the magnetic chemi cally 
peculiar stars in this tempe rature range (Rvabchikova 2005). We 
used the VALD database dKupka et alJI 19991) to compile a list 
of 8864 metal lines with the residual intensity exceeding 1% 
of the continuum in the 4000-6000 A wavelength region. The 
broad hydrogen Balmer lines and weak He i features were ex- 
cluded from the line list. The Zeeman splitting patterns required 
for the polarized spectrum synthesis and calculation of the LSD 
weights were obtained using Lande factors provided by VALD 
when available and estimated using the LS-coupling scheme oth- 
erwise. 

Theoretical spectra in four Stokes p arameters were calcu- 
lated on the adaptive wavelength grid (see lKochukhovl2007l) and 
then rebinned to an equidistant velocity spacing of 0.62 kms" 1 
and convolved with Gaussian profiles of FWHM - 2.5 kms 1 . 



1 The magnetic field vector is characterized by the field modulus B, 
angle y between the field vector and the line of sight and angle x, giving 
orientation of the field in the plane perpen dicular to the line of sight. 
Figure 2 in Piskunov & Kochukhov ( 2002) illustrates the definition of 
these angles. 
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Fig. 5. Small segment of the normalized local four Stokes parameter spectra calculated with synmast for the disk center (p. = 1), 
field strength B = 5000 G and magnetic vector orientation y = 45°, x - 0°- Linear and circular polarization profiles are shifted 
vertically for display purpose. Thin vertical bars plotted over the Stokes / spectrum indicate absorption lines deeper than 10% in the 
absence of magnetic field. 



This sampling and resolution corresp onds to the spec tra pro- 
vided by the HARPSpol instrument dSnik et all l2008h at the 
ESO 3.6-m telescope, which is the highest resolution stel- 
lar spectropolarimeter currently in operation. A few additional 
tests were done with the reduced resolving power of 65 000 
(FWHM -A.6 kms -1 ), corresponding to the polarization spec- 
tra from the ESPaDOnS and NARVAL instruments. 

Figure [5] shows a small segment of the local four Stokes pa- 
rameter profiles calculated with the help of synmast. 

3.2. Self-similarity of spectral lines 

First we examine the hypothesis of self-similarity of the Stokes 
profiles of different spectral lines, which is a fundamental as- 
sumption made by the LSD technique. We have selected all Fe i 
and Fe n lines from the original line list, removing a few iron 
lines with a very weak magnetic field sensitivity. The resulting 
1652 lines of diverse strengths and Zeeman patterns were re- 
distributed in wavelength to avoid blending. This line list was 
employed for the calculation of a grid of local Stokes parameter 
spectra for the field strength from 0. 1 to 5 kG and the field incli- 
nation from to 90°. The spectra were then convolved with the 
Gaussian instrumental profiles. 

A quantitative measure of the line self-similarity was ob- 
tained by computing the standard deviation of the line profiles 
interpolated onto a common velocity grid and divided by the 
appropriate LSD weights. We tested the standard LSD weights 
discussed in Sect. [2] as well as alternative weights obtained by 
substituting the central line depth d with the equivalent width 
expressed in velocity units, W v . 

The standard deviation of the rescaled Stokes profiles is plot- 
ted as a function of the field strength in Fig. [6] This figure shows 
results for two different resolutions. Clearly for Stokes / the line- 
self similarity assumption is a poor one and might be appropri- 
ate only for a comparatively low S /N of the observational data. 
Scaling by the equivalent width instead of the central line depth 
evidently gives better results. This reflects that strong lines are 
also systematically broader and therefore their central intensity 
is less affected by additional broadening. Based on this test we 



recommend to use the equivalent width rather than the line depth 
as an LSD weight in the applications concerned with Stokes / 
alone. However, the relative strength of the polarization signa- 
tures is described reasonably well by both d and W v . 



As is evident from Fig. [6] the assumption of the similarity 
of the line profile shapes is a reasonable one for Stokes V if the 
field strength does not exceed ~ 2 kG and the observational noise 
is > 5 x 10~ 3 (which corresponds to S /N - 200). Obviously, as 
the resolution and S /N increases, we see more and more differ- 
ence between individual spectral lines. The diversity of the linear 
polarization profiles is far greater than that of Stokes V. Even 
for the field strength of ~ 1 kG profiles noticeably differ from 
each other, and this difference rapidly increases with the field 
strength. Figure [6] shows that a definite reduction of the scatter 
between the scaled Stokes Q profiles is achieved by using G in- 

•^2 



stead of g in the corresponding line weight. 



The reason for the greater diversity of the Stokes Q line pro- 
files lies in their sensitivity to the line strength. As emphasized 
in Sect. 12.2.21 the linear polarization profiles are expected to 
be similar only when both the weak-field and the weak-line as- 
sumptions are satisfied, while for circular polarization only the 
former approximation is directly necessary. It is instructive to il- 
lustrate this by examining the average Stokes profiles for a very 
low field strength, when the weak-field approximation is valid. 
Figure Q shows that in this situation there is some difference be- 
tween the Stokes / line profile shapes of the weak and strong 
lines (as already discussed above) and a considerable difference 
between the shapes of Q profiles. The shape of the average linear 
polarization signature changes smoothly with the line strength, 
suggesting that the latter is the primary reason for the systematic 
discrepancy of the individual Stokes Q line profiles. In contrast, 
the Stokes V profiles are all very similar, almost independently 
of the line strength. Consequently, the basic assumptions of the 
LSD method are fairly robust for Stokes V but are more ques- 
tionable for linear polarization. 
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Fig. 6. Standard deviation of the scaled Stokes /, V and Q profiles as a function of the field strength for two different values of the 
field inclination y (45°, 90° for / and Q and 45°, 90° for V). Results for the two resolving powers are shown with red and blue 
curves. The solid, dashed, and dotted lines correspond to the standard deviations obtained using different line scaling factors. 
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Fig. 7. Mean scaled Stokes /, V, and Q signatures for the Fe lines of different central intensity. These profiles are obtained from the 
spectra computed with B = 100 G and an inclination angle y = 0° (Stokes / and V) and y = 90° (Stokes Q). 



3.3. The cutoff criterium and LSD profiles of individual 
elements 

In the second set of numerical experiments, conducted for the lo- 
cal non-magnetic spectrum, we investigated the effect of choos- 
ing a different cutoff criterium for the selection of lines that con- 
tribute to the LSD mask. Previous studies based on the LSD tech- 
nique arbitrarily applied a wide range of criteria: from a central 
depth of 0. 1 to 0.4. We found that for our noise-free spectra in- 
creasing the cutoff leads to the downward displacement of the 
resulting mean profile relative to the continuum, presumably be- 
cause of the combined effect of numerous weak lines that are 
excluded from the mask when the cutoff is too high. This con- 
tinuum depression of the LSD Stokes / profiles introduces an 
additional uncertainty in their interpretation. In subsequent tests 



we include in the mask all lines deeper than 10% in the absence 
of magnetic field, giving us a total of 3270 usable lines. 

We have also assessed the performance of our multiprofile 
extension of the LSD technique compared to the standard de- 
composition procedure. Using the non-magnetic spectrum we 
calculated the LSD profiles for Fe, Cr, and Ti in one case us- 
ing only the lines of each individual element and in another case 
accounting also for all other lines with the help of the second, 
"background" LSD profile. Figure [8] presents a comparison of 
the LSD profiles for the three elements obtained with these two 
varieties of the multiline technique. One can see that our multi- 
profile version of LSD is able to recover more reliable mean line 
profiles, free of the distortions caused by blends and without the 
spurious continuum offset. This improvement is most prominent 
for Ti, which has only 106 useful lines in our synthetic spectrum. 
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Fig. 8. Stokes / LSD profiles of individual elements recon- 
structed with the standard LSD technique (solid lines) and mul- 
tiprofile LSD (dashed lines). 

We therefore suggest the use of the multiprofile LSD in studies 
concerned with the analysis of small groups of lines present in 
stellar spectra. 

3.4. Longitudinal magnetic field 

In the next series of tests we studied how accurately the line of 
sight magnetic field component B z (longitudinal field) can be re- 
covered from the LSD Stokes V profiles. We computed a set of 
30 local polarized spectra, including spectral lines of all chem- 
ical elements, for a grid of magnetic field strengths (B = 100, 
200, 500, 1000, 2000, 5000 G) and the angles between the field 
vector and the line of sight (y = 0°, 22.5°, 45°, 67.5°, 90°). For 
each of these theoretical spectra we obtained the mean Stokes V 
profile and estimated the longitudinal field from its first moment, 
using Eq. (fJTJ. 

Figure [9] illustrates the relative error of B z as a function of 
the field strength and inclination. We find that the analysis of 
the LSD profiles recovers the true value of the longitudinal field 
with an accuracy of a few percent for the fields below 1 kG. At 
the same time, as is evident from Fig. [9] there is an uncertainty 
in the choice of the measurement window, resulting in another 
ss 2% error of B z . 

The longitudinal field estimate becomes less accurate when 
the field strength exceeds 1 kG. We find that for small y angles 
the LSD profile analysis overestimates the longitudinal field by 
up to 3%, while B z can be underestimated by as much as 10% 
for the intermediate magnetic field inclinations. 

We conclude that the derivation of the longitudinal magnetic 
field from LSD profiles appears to be quite robust but can be af- 
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Fig. 9. Relative error of the longitudinal magnetic field B z de- 
termined from the LSD Stokes V profiles. The lines of differ- 
ent colors show results for different inclinations of the magnetic 
field vector, while the solid and dashed lines illustrate the effect 
of choosing different integration limits. 



fected by the systematic errors of the order of a few percent. For 
strongly magnetic stars the longitudinal field tends to be under- 
estimated by the LSD technique. 

3.5. LSD profiles in Stokes parameters 

Another series of numerical experiments was carried out to an- 
swer the question whether or not the LSD Stokes parameter 
spectra behave as a real spectral line with respect to changes 
of the magnetic field strength and orientation. We derived LSD 
profiles from a grid of synthetic spectra covering the same pa- 
rameter space as studied in Sect. I3.4I The mask adopted in 
the least-squares deconvolution procedure was normalized us- 
ing (A) * Aq = 4700 A, d = (d) = 0.36 and g = (g) = 1.22. 
The resulting mean Stokes /, V and Q LSD profiles were com- 
pared with the theoretical spectra calculated for a fictitious Fe i 
line with the central wavelength A = 4700 A, an excitation po- 
tential of the lower level (3.47 eV) matching the average for this 
ion and a normal Zeeman triplet splitting pattern characterized 
by g = 1.22. The oscillator strength of the Fei line was adjusted 
to reproduce the equivalent width of the Stokes / LSD profile in 
the absence of magnetic field. 

The comparison of the LSD profiles and synthetic single-line 
calculations is presented in Fig. [I0]for a subset of the consid- 
ered magnetic field parameters. We found that even for the non- 
magnetic case the synthetic calculations for the resolving power 
R = 120000 (dotted lines in Fig. [TOjl do not match the LSD 
Stokes / line profile shape reconstructed from the full synthetic 
spectra having this resolution. The spectral resolution had to be 
decreased to approximately R = 80000 (solid lines in Fig. ITOb 
to reproduce the weak-field LSD profiles. This is probably an 
effect of the extended wings of strong lines, which contribute 
disproportionally to the mean profile. Thus, an estimate of the 
macroscopic velocity fields (stellar rotation, macroturbulent ve- 
locity) using LSD profiles might be affected by a small bias, es- 
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Fig. 10. Comparison of the LSD Stokes IQV profiles obtained from synthetic stellar spectra (symbols) with the theoretical line 
profiles computed for the average line parameters (lines). Different colors illustrate the results for three values of the magnetic field 
inclination with respect to the line of sight (0°, 45° and 90°). The four columns correspond to different magnetic field strength values 
(0.1, 0.2, 1 and 5 kG). The solid lines show theoretical profiles for the resolution R - 80000. The dotted lines in the Stokes / and 
V panels illustrate calculations fori? = 120000, which was the instrumental broadening applied to the original simulated stellar 
spectra. The dash-dotted lines in the Stokes V and Q panels represent calculations for the equivalent Zeeman splitting pattern with 
widely separated cr components. 



pecially for cooler stars in which the strong lines with devel- 
oped damping wings are considerably more important than for 
the T e ff = 9000 K stellar model considered here. 

Considering the behavior of the LSD Stokes I spectrum with 
the field strength, we found a fair agreement with the theoreti- 
cal single-line calculations up to B « 1 kG. Beyond this field 
strength value the lines start to exhibit partially resolved Zeeman 
patterns and the discrepancy between the LSD profile and its the- 
oretical single-line approximation becomes progressively worse. 

The circular polarization profiles show qualitatively the same 
picture as Stokes I. Below the field strength of ~ 1 kG there is an 
excellent agreement between theoretical Stokes V line shape and 
the Stokes V LSD profile, provided the extra Gaussian broad- 
ening is applied to the single-line calculations. The agreement 
becomes poor for B > 2 kG when individual Zeeman splitting 
patterns of spectral lines start to influence their Stokes V profiles. 
For strong fields, illustrated by the B = 5 kG profiles in Fig.fTUl 
the LSD Stokes V spectra are typically broader and have a lower 
amplitude compared to the single-line calculations. The discrep- 
ancy between LSD profiles and single-line calculations shows 
no clear dependence on the field orientation. 

In contrast to the encouraging weak-field behavior of the 
Stokes I and V, the Stokes Q single-line calculations fail in re- 
producing the corresponding LSD profiles for all field strengths. 
The shapes of the liner polarization signatures agree well for 



lower field strengths. But the amplitude of the LSD Stokes 
Q spectrum is systematically underestimated by the single-line 
synthesis for the fields weaker than ~ 1 kG. Then, as the field 
increases, the single-line Stokes Q amplitude grows faster than 
that of the LSD profile. The two amplitudes roughly agree for 
B * 1.5 kG and disagree in the other sense for stronger fields. 
This prominent field strength dependence implies that the dis- 
agreement that we revealed cannot be removed by a simple, 
field-strength-independent re-scaling of the LSD Stokes Q pro- 
files. 

The outcome of the calculations presented in this section 
suggests that in general the LSD profiles in four Stokes parame- 
ters do not behave as a single spectral line with a triplet Zeeman 
splitting pattern. It appears that the effective Lande factor and 
the line strength that successfully reproduce the LSD Stokes I 
and V spectra do not match the amplitude of LSD Stokes Q even 
for weak fields. There is also a significant systematic disagree- 
ment between the shapes of the LSD Stokes I and V profiles and 
single-line theoretical calculations for the field strengths above 
» 2 kG. 

In reality the majority of spectral lines exhibit anomalous 
Zeeman splitting, which can be considerably more complicated 
than a pure triplet. In principle, the positions and strengths of in- 
dividual components within the groups of red- and blue-shifted 
cr components can be manipulated to change the relative am- 
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plitude of the Stokes V and Q profiles, potentially providing a 
somewhat better agreement between the LSD polarization pro- 
files and single-line calculations. For example, the dash-dotted 
lines in Fig. [10] illustrate the Stokes Q profiles for a splitting 
pattern with widely separated pairs of cr+ and cr_ components. 
The effective Lande factor is kept at the value of 1 .22 and nei- 
ther the amplitude nor the shape of Stokes V change compared 
to the magnetic spectrum synthesis for the triplet splitting unless 
the field is very strong. On the other hand, the amplitude of the 
Stokes Q increases, providing a better match to the LSD Stokes 
Q profiles for a limited field strength range. However, parameters 
of this splitting pattern cannot be inferred from any properties of 
the LSD line mask and have to be calibrated empirically for a 
prescribed field strength range. 



3.6. Abundance determination from LSD profiles 

In this part of our investigation we assess to what extent LSD 
profiles behave as a real spectral line with respect to changes 
of abundances of individual elements. Specifically, we want to 
find out if there exists a set of atomic parameters reproducing 
changes of the LSD profile with abundance if it is modeled as 
a single spectral line. For the Ap-star like abundances adopted 
for our numerical tests the spectrum is dominated by the lines 
of Fe and Cr. Hence we considered the effects of the abundance 
variation for these two elements. 

Our LSD code allows to reconstruct an arbitrary number of 
independent LSD profiles using several masks. We took advan- 
tage of this feature in the abundance tests, reconstructing two 
LSD profiles: one for the element of interest and another for all 
other lines in the mask. 

Assuming the iron abundance [Fe]=+1 dex (hereafter abun- 
dances are given in the logarithmic scale with respect to the Sun), 
five non-magnetic theoretical spectra were calculated for vari- 
ous abundances of Cr in the range from +1 dex to +3 dex with 
a step of 0.5 dex. Then the chromium abundance was kept at 
[Cr]=+2 dex, while the abundance of Fe was changed from dex 
to +2 dex. This gave us five spectra with different Cr and five cal- 
culations with different Fe abundance, with a common spectrum 
for [Fe]=+1 dex and [Cr]=+2 dex. The latter spectrum was used 
for the selection of lines included in the LSD mask using our 
standard 10% cutoff criterium for the central line depth. 



3.6.1. Equivalent width 

The dependence of the equivalent width (EW) of the Fe LSD 
profile on the abundance of this element is illustrated in Fig.fTTI 
We tried reproducing this dependence with single-line calcula- 
tions using different methods of determining the parameters of 
the average line. Both the Fei and Fen lines contribute at a 
similar level to our synthetic spectra and to the resulting LSD 
profiles. Hence one can adopt mean line parameters by averag- 
ing over the lines of particular ion or averaging over both ions. 
Furthermore, one can determine the mean oscillator strength us- 
ing the harmonic or arithmetic average, the latter giving system- 
atically higher value of the oscillator strength. This yields six 
possible combinations of ion, excitation potential, and oscillator 
strength. It is not clear which one should be chosen for modeling 
LSD profiles, hence we tested all of them. The resulting equiva- 
lent width changes with abundance are illustrated with the lines 
in Fig.QT] Evidently, two or three combinations of the mean line 
parameters give somewhat better results, but all predict much 



0.5 

? °-° 
--o, 

w> -1.0 
-1.5 
-2.0 




Fe II (E=7.59 eV < 
Fe II (E,=7.59 eV < 



log gf= 
log gf= 



0.0 



0.5 



1.0 
[Fe] (dex) 



1.5 



2.0 



Fig. 11. Dependence of the equivalent width on abundance for 
the LSD profiles (symbols) and theoretical profiles of Fe spec- 
tral lines (lines). Different lines illustrate the EW-abundance de- 
pendence obtained for different combinations of line parameters. 
Red: averaging over all Fe lines; green: averaging over Fe i lines; 
blue: averaging over Fe n lines. The dashed and solid lines show 
results for the oscillator strength determined using harmonic and 
arithmetic averaging, respectively. 
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Fig. 12. Comparison of the Fe LSD profiles (symbols) and the 
profiles of theoretical spectral line (solid curve) for one of the 
sets of atomic parameters used in Fig.fTTI 



steeper slopes of the EW dependence on abundance than shown 
by the LSD profile. 

Calculations for Cr show a very similar discrepancy and re- 
veal the same pair of the line parameters producing a slightly 
better fit to the LSD EWs compared to other options. 

3.6.2. Line profiles 

We complement the EW comparison with a more detailed in- 
vestigation of the Fe LSD profiles as a function of abundance. 
This could potentially shed some light on the discrepancy be- 
tween the LSD and expected single-line dependence on element 
abundance. 

As shown in Fig. [12] the LSD profiles appear to be stronger 
than the synthetic spectrum for lower abundances and weaker, 
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but wider, for larger abundances. These broad wings do not con- 
tribute significantly to EW, and their exclusion would only de- 
grade even more the agreement between the LSD profile and 
single-line abundance behavior. 

We have verified that blending of the Fe or Cr lines by the 
lines of other elements is not the reason for the observed abun- 
dance dependence. In particular, if we calculate the pure Fe or 
Cr line spectrum and de-blend all lines by shifting them in wave- 
length, we still get the same slope of the EW vs. abundance trend 
as illustrated in Fig.fTTl 

We conclude that weak and strong lines are combined by the 
LSD technique in a way that the resulting mean Stokes / profile 
does not show the abundance dependence of a real spectral line. 
This is probably because of a highly non-linear change of both 
the central depth and the width of stronger absorption lines with 
abundance variation. This suggests that deducing the absolute 
elemental abundances from LSD profiles is a highly non-trivial 
task, which can be accomplished only using LSD profile calcu- 
lations from synthetic spectra similar to those presented here. 
Nevertheless, LSD profiles can still be used for a comparison of 
the element concentrations between different stars with similar 
fundamental parameters or different regions on the surface of the 
same star. 

4. Conclusions 

We have carried out a detailed study of the least-squares decon- 
volution of both intensity and polarization stellar spectra. This 
multiline signal-enhancing technique is widely used in a broad 
range of modern spectroscopic and spectropolarimetric studies, 
but lacks a thorough justification and was never tested using 
synthetic spectra. Our study fills this gap, clarifying the basic 
assumptions of LSD, investigating various ways to improve its 
numerical implementation and, finally, testing the sanity of the 
physical information derived from LSD spectra. 

The main conclusion of our study can be summarized as fol- 
lows. 

1 . The basic LSD assumptions of line profile self-similarity and 
linear addition of blends are satisfied only for very weak 
(less than 20% of the continuum) absorption lines in Stokes 
/. Similar results are obtained for linear polarization signa- 
tures, which exhibit a prominent line-strength dependence in 
weak fields and become quite diverse in shape owing to dif- 
ferent Zeeman splitting patterns for ^1 kG magnetic fields. 
On the other hand, the LSD assumptions are appropriate for 
the majority of lines in Stokes V for weak and moderately 
strong magnetic fields (B < 2 kG). 

2. We developed a new, improved LSD code, which is general- 
ized to handle multiple mean profile reconstruction from the 
same spectrum and employs regularization to suppress ex- 
cessive numerical noise that appears in the LSD profiles as 
a result of the autocorrelation matrix inversion. At the same 
time, we demonstrated th at the recently suggested non-linear 
deconvolution method of Sennhaus er et al.l (12009) is applica- 
ble only to the spectra of very slowly rotating stars observed 
at the highest available spectroscopic resolutions, but is in- 
ferior to the standard LSD method when the line blending is 
primarily due to rotational and instrumental broadening. 

3. We emphasized the importance of choosing correct line 
weights in the calculation of LSD profiles. It is the normal- 
ization of these weights, but not the mean line parameters, 
that determine the characteristics (central wavelength, Lande 
factor, etc.) of the resulting LSD profile. We suggested an 



improved weighting scheme for the computation of the linear 
polarization LSD profiles, which accounts for the anomalous 
Zeeman splitting of spectral lines. 

4. We performed comprehensive numerical experiments in 
which we analyzed LSD profiles computed from the local 
theoretical Stokes parameter spectra of a mid-A magnetic 
chemically peculiar star. The measurements of the mean 
longitudinal magnetic field from the synthetic LSD profiles 
show that B z accurate to within a few % can be retrieved for 
the magnetic field strength of up to ~1 kG, but systematic 
errors of up to 10% can be expected for stronger fields. 

5. Generally, the LSD intensity and polarization profiles do not 
behave as a spectral line with parameters determined by av- 
eraging over individual lines contributing to the mean spec- 
trum. In particular, the LSD Stokes / profiles appears broader 
and its equivalent width changes to a lesser extent with vari- 
ation of element abundances than for a real spectral line. The 
LSD Stokes V and Q cannot be described simultaneously by 
the same line parameters even in a weak magnetic field. The 
only parameter region where we found the single-line ap- 
proximation of the LSD profiles to be reasonably robust is 
the circular polarization spectra in the 0-2 kG field strength 
interval. Whenever the field intensity exceeds this limit or 
an element abundance is changed or Stokes Q is considered 
simultaneously with the circular polarization, the LSD pro- 
file cannot be predicted by a single-line spectrum synthesis. 
These results suggest that the single-line approximation of 
LSD profiles should be abandoned in favor of a direct com- 
parison of the mean profiles inferred from observations and 
from the same lines in theoretical polarized synthetic spec- 
trum. 
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